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Abstract 

This is a theoretical study of the interplay of optical phase-conjugation and 
multiple scattering. We calculate the intensity of light reflected by a phase- 
conjugating mirror when it is placed behind a disordered medium. We com- 
pare the results of a fully phase-coherent theory with those from the theory 
of radiative transfer. Both methods are equivalent if the dwell time Tdweii 
of a photon in the disordered medium is much larger than the inverse of 
the frequency shift 2Aa; acquired at the phase-conjugating mirror. When 
TdwcuAo; < 1, in contrast, phase coherence drastically affects the reflected 
intensity. In particular, a minimum in the dependence of the reflectance on 
the disorder strength disappears when Aw is reduced below l/rdweii- The 
analogies and differences with Andreev reflection of electrons at the interface 
between a normal metal and a superconductor are discussed. 

PACS numbers: 42.65.Hw, 42.68.Ay, 42.25.Bs, 78.20. Ci 
physics/9708010 
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I. INTRODUCTION 



Phase conjugation is the reversal of the sign of the phase of a wave function. A phase- 
conjugated wave retraces the path of the original wave, thereby canceling all accumulated 
phase shifts. Phase conjugation was first discovered for electronic waves [|l|, and later for 
optical waves [0,0. For electrons, phase conjugation takes place at the interface between 
a normal metal and a superconductor. An electron at energy E above the Fermi energy 
Ep is reflected at the angle of incidence (retro-reflected) as a hole at energy E below Ep, 
a process known as Andreev reflection A phase-conjugating mirror for light consists of 
a cell containing a liquid or crystal with a large nonlinear susceptibility, pumped by two 
counter-propagating beams at frequency ujq. A wave incident at frequency c^o + Acj is then 
retro-refiected at frequency uq — Auj, a process known as four- wave mixing 

The interplay of multiple scattering by disorder and phase conjugation has been studied 
extensively in the electronic case, both experimentally and theoretically. (See Ref. for 
a review.) In the optical case the emphasis has been on weakly disordered media, which 
do not strongly scatter the waves 0. Complete wave- front reconstruction is possible only 
if the distorted wave front remains approximately planar, since perfect time reversal upon 
reflection holds only in a narrow range of angles of incidence for realistic systems. (For the 
hypothetical case of perfect time- reversal at all angles, see Ref. |T^.) McMichael, Ewbank, 
and Vachss measured the intensity of the reconstructed wave front for a strongly inho- 
mogeneous medium (small transmission probability Tq), and found that it was proportional 
to Tq — in agreement with the theoretical prediction of Gu and Yeh If Tq <C 1, the 
intensity of the reconstructed wave is much smaller than the total reflected intensity. The 
total reflected intensity was not studied previously, perhaps because it was believed that the 
diffusive illumination resulting from a strongly inhomogeneous medium would render the 
effect of phase conjugation insignificant. In this paper we show that a strongly disordered 
medium backed by a phase- conjugating mirror has unusual optical properties, different both 
from the weakly disordered case and from the electronic analogue. 

We distinguish two regimes, depending on the relative magnitude of the frequency shift 
2Au! acquired at the phase-conjugating mirror and the inverse of the dwell time Tdweii of 
a photon in the disordered medium. (For a medium of length L and mean free path /, 
with light velocity c, one has Tdweii — L? /cl.) In the coherent regime, Acj ^ 1/Tdweii, phase 
conjugation leads to a constructive interference of multiply scattered hght in the disordered 
medium. In the incoherent regime, Aa; ^ 1/rdweii, interference effects are insignificant. In 
both regimes we compute the reflectances i?+ and defined as the reflected power at 
frequency ujq ± Au; divided by the incident power at frequency ujq + Acj. A distinguishing 
feature of the two regimes is that (in a certain parameter range) the reflectance i?_ decreases 
monotonically as a function of L/l in the coherent regime, while in the incoherent regime it 
first decreases and then increases. 

The outline of this paper is as follows. After having formulated the problem in Sec. 0, we 
discuss in Sec. |T| its solution using the Boltzmann equation, ignoring phase coherence. This 
is the theory of radiative transfer |T3|JT^ . A simple result is obtained if we neglect angular 
correlations between the scattering in the disordered medium and at the phase-conjugating 
mirror. We compare this approximation with an exact solution of the Boltzmann equation. 
In Sec. 1^ the phase-coherent problem is addressed, analytically using random-matrix theory, 
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FIG. 1. Schematic drawing of the disordered medium backed by a phase-conjugating mirror 
(PCM). Light incident at frequency loq + Alo is reflected at the two frequencies loq ± Auj. 



and numerically using the method of recursive Green functions. Results of this section were 
briefly presented in Ref. |]T^. We conclude in Sec. |V| with a comparison with the electronic 
analogue of this problem. 



II. FORMULATION OF THE PROBLEM 

We study the system shown in Fig. |l]. It consists of a disordered medium (length L, mean 
free path /), backed at one end by a phase-conjugating mirror. The other end is illuminated 
diffusively at frequency Uj^ = uq + Au, where ujq is the pump frequency of the mirror. We 
are interested in the amount of hght reflected at frequency and = Uq — Au. 

To reduce the problem to the scattering of a scalar wave, we choose a two-dimensional 
geometry. The scatterers consist of dielectric rods in the z-direction, randomly placed in 
the x-y-p\a.ne. The electric field points in the z-direction and varies in the x-y-plane only. 
Two-dimensional scatterers are somewhat artificial, but can be realized experimentally |]16| . 



We believe that our results apply qualitatively to a three-dimensional geometry as well, 
because the randomization of the polarization by the disorder renders the vector character 
of the light insignificant. 

The z-component of the electric field at the frequencies Uj^ and U- is given by 

E±{x, y, t) = Re S±{x, y) exp(— ic<j±t). (2.1) 

The phase- conjugating mirror (at x = L) couples the two frequencies via the wave equation 

mm 



Ho 7* \ fS+\ 2eAuj f S- 



7 -nj\£*_j coo \£ 



(2.2) 



The complex dimensionless coupling constant 7 is zero for x < and for x > Lc, with Lc 
the length of the nonlinear medium forming the phase-conjugating mirror. For < x < Lc 
it is proportional to the electric fields £1, £2 of the two pump beams and to the third-order 
nonlinear susceptibility Xs- 

7 = -7^X3^1*^2 = loe'"^, 0<x<Lc. (2.3) 

ZEo 



The Helmholtz operator Ho at frequency uo is given by 



Ho = -fco"'V^ - £, (2.4) 

where e{x,y) is the relative dielectric constant of the medium. We take e = 1 except in the 
disordered region —L < x < 0, where e = l + 6e{x,y). The fluctuations 6e lead to scattering 
with mean free path /. We assume k^l ^ 1, where ko = uJo/cis the wave number of the light 
(velocity c). The validity of Eq. ( p^ ) requires Au/uo <^ 1 and I7I =70 <^ 1- The ratio of 
these two small parameters 

. 2Aa; , , 

6 = (2.5 

is a measure of the degeneracy of the incident and the reflected wave, and can be chosen 
freely. 

In the absence of disorder, an incoming plane wave in the direction (cos 0, sin 0) is retro- 
reflected in the direction (— cos0, — sin0), with a different frequency and amplitude. The 
scattering matrix for retro-refiection is given by ||5|JT7|-pO 



£+y_f -m(0)e-*'^\ /£ 



£1J \ia{(p)e'^ J \£ 



[2.6a) 



a(0) = [Vl + 52cotan {aV 1 + 6^ / cos (f)) + i6]-\ (2.6b) 
a = l^okoLc. (2.6c) 



The crucial difference with Ref. [0] is that the reflectance is angle dependent and that the 



reflection matrix is non-Hermitian. This implies that not all phases will be canceled in the 
conjugation process. In Fig. we have plotted the reflectance | function of the angle 

of incidence for a = ir/A and two values of 5 = 0.75 and 0.9. The value a = 7r/4 is chosen 
such that a = 1 for normal incidence at frequency ujq (i.e. for = 0,5 = 0). The two values 
of 6 have been chosen such that the angular average of the reflectance, 

A = r'"^ d(t) COS0 |a(0)|^ (2.7) 



is > 1 for 5 = 0.75 and < 1 for 6 = 0.9. (The cos0 weight factor in Eq. (|2.7| ) corresponds 
to diffusive illumination.) In most of the numerical examples throughout this paper we will 
use these values of a and 6. 



III. PHASE-INCOHERENT SOLUTION 

A. Radiative transfer theory 

Within the framework of radiative transfer theory |T^|1^, the stationary distribution 
I{x,y,(j)) oc \£\'^ of the light intensity, at frequency u and wavevector (A;cos0, fcsin0), is 
governed by the Boltzmann equation 

( d d\ 1 /•2'r 

Mcos0— + /sin0— 1 /(x,|/,0) = -/(a;,?/,0) + — d(j)' I{x,y, (f)'). (3.1) 
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FIG. 2. Reflectance of the phase-conjugating mirror as function of the angle of incidence, 
computed from Eq. ( |2.6D for two choices of parameters. 




We neglect absorption and assume isotropic scattering in the x-y plane, with mean free 
path /. The phase-conjugating mirror couples the intensities J± of light at the two frequencies 
ci;± = cuo ± Acu. We assume that / is independent of frequency. The symmetry of the system 
implies that J(x, 0) = /(x, |0|). In this section we take e [0, vr]. For each frequency the 
Boltzmann equation takes the form 



/ cos I 



dx 

I±{x) 



J±(x) 
1 

TC Jo 



I±{x, 

')I±{X, 



(3.2a) 
(3.2b) 



Eq. ( p.2| ) has to be supplemented by boundary conditions at the two ends x = —L and 
X = of the disordered medium. We consider a situation that the system is illuminated at 
X = —L with diffusive light at frequency uj+, hence 



I+{—L,(j)) = Iq, for COS0 > 0, 
J-f-L, (/)) = 0, forcos(/)>0. 



(3.3a) 
(3.3b) 



At a; = the light is reflected by the phase-conjugating mirror. The intensity is multiplied 
by 



|a(0)p 



sin^(aVl + 5^/ cos0) 
5^ + cos^(aVT+5^/ cos< 



(3.4) 



according to Eq. (|2.6|) . The reflection is accompanied by a change in frequency uj± 
so that the boundary condition is 

/±(O,0) = |a(0)p/^(O,7r - 0), for cos0 < 0. 

The flux j± associated with the intensity I± is defined by 

/■TT 

j± = / d(j) COS0 J±(a;,0), 



UJ: 



(3.5) 



(3.6) 



and is independent of x [dj±/dx = according to Eq. ( p.2|) ]. The reflectance R- is defined 
as the ratio of the outgoing flux at frequency c<j_ and the incoming flux at frequency uj+, 
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i?_ = -j_/Jo. (3.7) 
The total outgoing flux is (i?_ + R+)Io, where 

R+ = l-u/Io (3.8) 
is the ratio of the outgoing flux and the incoming flux at the same frequency u^. 

B. Neglect of angular correlations 

A simple analytical treatment is possible if the angular correlations between multiple 
reflections by the disorder and the phase-conjugating mirror are neglected. Here we present 
this simplified treatment, and in the next subsection we compare with an exact numerical 
solution of the Boltzmann equation. 

We first consider the disordered region by itself. The plane-wave transmission probability 
|t((/))P is the ratio of transmitted to incident flux when the incident light is a plane wave in 
the direction (cos 0, sin 0). The transmission probability T for diffusive illumination is then 
given by 

tt/2 

ci0 cos0|^(0)|^ (3.9) 



such that T is the fraction of the flux incident from a diffusive source which is transmitted 
through the disordered region. This probability has been calculated in Ref. PT| from the 
Boltzmann equation (|3.2|) . The result is 

T = {l + 2riL/7iiy\ (3.10) 

where t] is a numerical coefficient which depends weakly on L/l. In the ballistic limit 
{L/l 0) 7] has the value vr^/S and in the diffusive limit {L/l oo) r] equals 1. In this 
subsection (but not in the next) we take rj = 1 for all L/l for simplicity. 

We use Eq. ( p.lOp to obtain the reflectance R± for the case that the disordered medium 
is backed by a phase-conjugating mirror with reflectance 

A = dc^ cos0sin^(a.^rTF/cos0) ^^^^^^ 
Jo 52 _^ pQg2 {ayl + 5"^/ cos 0) 

Since T and A are angular averages, we are neglecting angular correlations. The light 
that comes out at frequency uj- has been reflected an odd number of times at the mirror. 
The light that has been reflected once has traversed the medium twice, which leads to a 
contribution T'^A to -R_. Light that has been reflected three times by the mirror contributes 
T^yl'^(l— T)^, since it has been reflected two times by the medium (each time with probability 
1 — T) . Summing all contributions, one finds 

R. = T^A + T^A^l - T)2 + T^A'il - T)^ + ■ ■ ■ = /^'^ (3.12a) 

1 — (1 — i j A 
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Light that comes out at frequency ui^ has been reflected an even number of times at the 
mirror. Zero reflections by the mirror contributes 1 — T to two reflections contributes 
T^yl^(l — T), and four reflections T^A^(1 — T)^. Summing the series, one finds 



T + 



T2(l -T)A2 



(3.12b) 



The geometric series leading to Eq. (p.l2|) diverges if (1 — T) A > 1. This indicates that 
there is only a stationary solution to the Boltzmann equation if both the gain at the mirror 
and the scattering in the medium are sufficiently weak. If A is increased at fixed a = 7r/4 
by reducing 6, the refiectances R± diverge when 6 = 6c- (This divergence is preempted by 
depletion of the pump beams in the phase-conjugating mirror, which we do not describe.) In 
the approximation of this subsection, Sc is determined by {1—T)A = 1, or L/l = ^n{A — l)~^. 
In the ballistic limit, T = 1 and A < oo for any 6 > 0. In the diffusive limit, T = and 
A = 1 foT 6 = 0.78. Hence, 6c increases from to 0.78 as L/l increases from to oo. 



C. Exact solution of the Boltzmann equation 



The Boltzmann equation (|3.2|) can be solved exactly, without neglect of angular correla- 



tions, by adapting the method of Ref. to an angle-dependent boundary condition. We 
first rewrite Eq. ( p.2|) as 



1 



|-e^/'™^^J±(x,0)= , 
ox I cos 



(3.13) 



and then integrate once over x, using the boundary conditions ( |3.3| ) and ( |3.5|) . The result 
is 



/+(x, 



I-(x, 



I-(x, 



^ dx' 
-L I COS <; 
dx' 



I COS( 



) 



+ e"^/'™" "^10(0)12 /_(0,7r- 



^ dx' 
-L I COS <; 
dx' 



I COS( 



. ^-(x-x')/lcos (t> J ^^/^ 
) 

^-{x-x')/lcos<t> j^l^^/-^ 



+ e-"/'"°^<^|a(0)|2/+(O,7r- 



for COS0 < 0, 
for cos (p > 0, 

for cosd) < 0. 



(3.14a) 

(3.14b) 
(3.14c) 

(3.14d) 



Substitution of Eqs. ( p.l4c|) and (|3.14a| ) into, respectively, Eqs. ( |3.14b|) and ( |3.14d| ) yields 



I+(x, 



dx' 



X Z|cOS( 



X r -i^e^'/'l™^*l/_(x') 

L l\ COS© 



for coscj) < 0, (3.15a) 
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I_{x,(f))= / „ , e-("'-")/'l^°'^'^l I_{x') + e"/'l^°^'^l|a(0)|' 

Jx l\cOS(f)\ 

j^g-L/«[ COS 01 ^ / ^ -e^'/''=°^'^l/+(x') , forcos0<O. (3.15b) 

J~L l\ COS 01 J 

Finally, integration over cj) leads to two coupled integral equations for the average intensities, 



lAx 



1-0 - /"O 

j ^dx' Mi{x,x') I+{x') + J ^dx M2{x,x')I^{x) + Qi{x)Io , (3.16a) 

I_{x)= r dx' Mi{x,x')I_{x')+ r dx' M2{x,x')I+{x') + Q2{x)Io. (3.16b) 

J-L J-L 



We have defined the following kernels and source terms: 

Mi(x,x') = - r^'^l^e-l^-^'l/'™^'^ = liTodx-x'l/O, (3.17a) 

TT JO / COS TT/ 

M2(x,x') = - Jt-e^-+-')n-os<t> |^(^)|2^ (3 17b) 

TT JO f COS (j) 

Q^{x) = 3 /"^^ci0e-(^+^)/'^°^<^, (3.17c) 

-{L-x)/lcos^ I / ,\|2 



TT Jo 

1 /-Va 



g2(a;) = - / rf0e-^^-"^/'^°^<^ |a(0)|^ (3.17d) 

TT JO 

where is a Bessel function. 

Eq. (|3.16|) is the analogue for the present problem involving two coupled frequencies 
of the Schwarzschild-Milne equation in the theory of radiative transfer |]T3|,|H|. We have 
solved it numerically by discretizing with respect to x so that the integral equation becomes 
a matrix equation. From the average intensities I±{x) one finds the intensities I±{x,(f)) 
using Eqs. ( |3.14| ) and (|3.15|) . The refiectances R± then follow from Eqs. ( p.6| )-( p^ . For 



numerical stability we have imposed a cut-off on the rapidly oscillating function a(0) at 
grazing incidence, by setting a(0) = for 0.497 vr < < |7r. 

In Figs. ^ and ^ we show results for I±{x) and R± for a = tt/A and 6 = 0.75 and 0.9. For 
S = 0.75 there is an effective gain at the mirror {A > 1), while for 6 = 0.9 there is an effective 
loss {A < 1). For an ordinary mirror one can show that /±(0) = |Jo. Instead, we find that 
/_(0) > /+(0) > iJo for 6 = 0.75, indicating gain, and J_(0) < /+(0) < ^Iq for 5 = 0.9, 
indicating loss. In each case the density profiles are approximately linear in the bulk, with 
some bending near the boundaries at a; = —L and x = 0. For 6 = 0.75, both /2_ and -R+ 
diverge when L/l = 28, while for 6 = 0.9 no such divergence occurs. As discussed earlier, 
the divergence indicates that for 6 = 0.75 and L/l > 28 there is no stationary solution to 
the Boltzmann equation. For fixed L/l and a, the divergence of R± occurs at a critical 
value 6c, such that a stationary solution requires 6 > 6c. The dependence of 6c on L/l at 
fixed a = 7r/4 is plotted in Fig. ||. 

In Figs. § and |^ we also compare the exact numerical solution of the Boltzmann equation 
of this subsection with the approximate analytical solution ( ^.12| ) of the previous subsection. 
As one can see, the agreement with the exact results is quite good. 
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FIG. 3. Intensity profiles in the disordered medium, computed from the exact numerical solu- 
tion of the Boltzmann equation, for a = 7r/4 and two values of 6. (a) is for a nearly ballistic system 
{L/l = 1), (b) is for a diffusive system {L/l = 15). 
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L/l 

FIG. 4. Reflectance i?_ as function L/l, computed from the exact solution of the Boltzmann 
equation for a = n/A and 5 = 0.75 (dashed curve), 5 = 0.90 (solid curve). The dotted curves are 
the approximate result ( p.l2a| ), in which angular correlations are neglected. The inset shows the 
exact reflectances R± for 6 = 0.75, over a broader range of L/l (logarithmic scale). For 5 = 0.75 
the reflectances diverge at L/l = 28. No divergence occurs for 6 = 0.90. 
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FIG. 5. A stationary solution to the Boltzmann equation requires 6 > 5c- The solid curve is the 
exact result for 5c (at fixed a = vr/4, as function of L/l), the dotted curve follows from Eq. ( p. 12 ), 
obtained by neglecting angular correlations. 
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IV. PHASE-COHERENT SOLUTION 



A. Scattering matrices 

We now turn to a phase- coherent description of the scattering problem. To define finite- 
dimensional scattering matrices we embed the disordered medium in a waveguide (width W), 
containing N± = liat{uj±W/c'7i) 3> 1 propagating modes at frequency A basis of scatter- 
ing states consists of the complex fields 

E^ „{x, y, t) = k^]l^ sin (^^) exp(ifc±,„x - iuj±t), (4.1a) 

E±_„(x,y,t) = k^]l^ sin (^^) exp(-iA;±,„x - iuj±t). (4.1b) 

Here n = 1, 2, . . . , N± is the mode index and the superscript > (<) indicates a wave moving 
to the right (left), with frequency uj± = ujq ± Acj and wavenumber 

k±,n = (^l/c' - n'TT^/W^f/^ (4.2) 

The normalization in Eq. ( [4.1|) has been chosen such that each wave carries the same fiux. 
With respect to the basis ( [4.1|) , incoming and outgoing waves are decomposed as 



E^'^ = U+,nEX^n + E ^-,nE^,n. (4.3a) 
n=l 71=1 

N+ Af_ 

= E ^+,nEt^n + E V-,nE<^^. (4.3b) 
n=l n=l 

The complex coefficients are combined into two vectors 

u = (-u+^i, • • • , 'u+,7v+, _i, u*_ 2, . . . , Ml,iv_)'^, (4.4a) 

V = (w+,i,w+,2, • • . ,w+,Ar+,w*_i,t;*^25 • • ■,V-,nJ^- (4-4b) 
The refiection matrix r relates u to v, 

v = ru, r ^+-). (4.5) 

The dimension of r is (A^+ + iV_) x {N^ + A^_), the submatrices r±^-|- have dimensions 
N± X A^-j-. For Aa; ^ cuq we may neglect the difference between iV_|_ and N_ and replace 
both by AT = lnt(A;ol^/7r). 

In the absence of disorder the refiection matrix is entirely determined by the phase- 
conjugating mirror, 

/ -tae-'^ \ 

rpcM=(^^^^.^ j' ^^-^^^ 

amn = a{(f)n)Smn, 0„ = arcsin(n7r/ /cq W^) ■ (4.6b) 



The elements of the N x N diagonal matrix a are obtained from Eq. ( |2.6| ) upon substitu- 
tion of (j) by 0„, being the angle of incidence associated with mode n. (The difference in 



11 



angle between the two frequencies ojj^ and U- can be neglected if Acu <^ ujq.) The angular 
average ( p.7|) of the reflectance corresponds to the modal average 

A = — Tr aa^ (4.7) 

In the limit N ^ oo the two averages are identical. 

The disordered medium in front of the phase-conjugating mirror does not couple a;+ 
and Its scattering properties at frequency u are described by two N x N transmission 
matrices t2i{uj) and ti2(ix;) (transmission from left to right and from right to left) plus two 
N X N reflection matrices rii(u;) and r22{uj) (reflection from left to left and from right to 
right). Taken together, these four matrices constitute a 2N x 2N scattering matrix 

Sdisordcr(^)-(^^^^(^) r22(c.)j' ^^-^^ 

which is unitary (because of flux conservation) and symmetric (because of time-reversal 
invariance). It is simple algebra to express the scattering matrix r of the entire system in 
terms of the scattering matrices rpcM and Sdisordcr of the phase-conjugating mirror and the 
disordered region separately. The result is 

r++ = r-^i{uj^) + ti2(w+)ar;2(^_)a[l - r22(a;+)ar;2(^_)a]~H2i(a;+), (4.9a) 

r__ = rl^{u_) + ti2(^_)ar22(^+)a[l - r;2(^_)ar22(^+)a]"H;i(^_), (4.9b) 

r_+ = ie'^tl2{uj_)a[l - r22{^+)air;2{uj_)a]-\,{uj^), (4.9c) 

r+_ = -ze-'%2i^+W - r;2(^-)ar22(^+)a]-H;i(cu_). (4.9d) 

We seek the reflectances 

R- = ^Trv_^rU; = ^Tr r++rt+, (4.10) 



averaged over the disorder. We will do this analytically, using random-matrix theory 
and numerically, using the recursive Green function technique |2^. We consider two differ- 
ent regimes, depending on the relative magnitude of Aa; and 1/rdweib where Tdweii — L'^/d is 
the mean dwell time of a photon in the disordered medium. If TdweiiAu; <^ 1 the difference 
between Sdisorder(i^+) and Sdisorder('^-) is insignificant, because the phase shifts accumulated 
in a time Tdweii are approximately the same for frequencies Uj^ and u^. We call this the 
coherent regime. If Tdweii Acu ^ 1, on the contrary, phase shifts at and are essen- 
tially uncorrelated, so that Sdisorder('^+) and Sdisordcr (t^-) are independent. We call this the 
incoherent regime. 

B. Random-matrix theory 

Without loss of generality the reflection and transmission matrices of the disordered 
region can be decomposed as [E^ 



r,,{uj^) =iY^^l - Vl, r22(a;±) =iV^^jl - T± U^, (4.11a) 
t,2(^±) = V± Ul, t2i(^±) = \J^JtI Yl (4.11b) 
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Here U± and V-t are NxN unitary matrices (we take = = in this subsection) and 
T± is a diagonal matrix with the transmission eigenvalues r± „ G [0, 1] on the diagonal. The 
subscript ± refers to the two frequencies u!+ and . In this so-called "polar decomposition" 
the reflectances R± take the form 

1 



= ^Tr T_0 f 1 - a/1 - T+ -T ft 



■ T+ fl - n^Jl-T_n*Jl-T+] ^ n\ (4.12a) 



R+ = ^Tr (1 - Th 



-1 



-1 



• T , ( 1 - ^i^/l - T_ n*Ji - T , nU/1 - t_ fi 



1 



- ^Tr T^ ji - ( 1 - o^/l - T_ n\/i - 1 ^^^/l - T_ n* 



-1 



1 



-1 



- — Tr T+^1 -T+f^Yl - T-f^ (^1 - - T+f^ ~ '^-^j ' (4.12b) 
n = UlaU+. (4.12c) 

To compute the averages {R±) analytically in the large A^-limit we make the isotropy 
approximation ^2 that the matrices U± and V-t are uniformly distributed over the unitary 
group U{N). This approximation corresponds to the neglect of angular correlations in 
the radiative-transfer theory (Sec. PIB|) . For Xdweii'^'^ ^ 1 we may identify U+ = U_ and 
V+ = V_. For TdweiiAcoi ^ 1 we may assume that U+, U_, V+, and V_ are all independent. 
In each case the integration / (iU/(U) over U{N) with ^ 1 can be done using the large 
iV-expansion of Ref. [El]. The remaining average over r-|-„ can be done using the known 



density p(r) of the transmission eigenvalues in a disordered medium p2 



The calculation is easiest in the incoherent regime (rdwcuAct; ^1). The integration over 



IA{N) can be carried out using the formula ||24 



jd\J j rfV^Tr (A1UA2VA3U ■ • • Ap)(BiUB2VB3U ■ • • B,)t = 5„N-p f{ Tr A.Bf 

+ 0{N-P-^). (4.13) 

To apply this formula we expand the inverse matrices in Eq. ( [4.12| ) in a power series in U± 
and integrate term by term over the independent matrices U_|_ and U_. The result is, to 
leading order in A^, 

/ rfU_ / rfU+ R- = Y. T-T+A^P+^l - T_)P(1 - T+Y 

~ 1- (1-T_)(1-T+)A2' ^^-^^"^^ 

. i?+ = 1 - T+ + — , (4.14b) 

+ + + 1 - (1 -T_)(l -T+)A2' ^ ' 
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where we have defined the modal average 



1 1 ^ 

T± = j;^TrT^ = -Y.J±,n. (4.15) 



The modal average A was defined in Eq. ( [4.7|) . The quantities T± still depend on the 
configuration of the scatterers, but the fiuctuations around the average {T±) are smaller by 
an order 1/A^ than the average itself. Moreover, the average (T±) equals the transmission 
probability T from radiative transfer theory, Eq. ( |3.10| ), again up to corrections of order 
1/A^. Replacing T± in Eq. ( [4.14| ) by T we obtain 

(i?.) = l-T + ^-i^-^, (4.16b) 



which is the result (p. 12 ) of radiative transfer theory with neglect of angular correlations. 



The conclusion is that in the incoherent regime phase coherence has no effect on the re- 
flectance of the system to leading order in A^. 

The situation is entirely different in the coherent regime (rdwei!'^'^ -C 1). To see the 
difference it is instructive to first consider the simplified model that the matrix amn = CLo^mn 
is proportional to the unit matrix (a scalar). Because U_ = U+ for Tdweii^i^ 1, we then 
have fimn = cio^mn- There is therefore no average over U{N) to perform. We only have to 
average over one set of transmission eigenvalues r+_„ = r_^„ = r„. This average amounts to 
the integrals 



(fl-> = i|*pM^J|tl^. (4,17a) 



^ -'o |i _ _|_ 

The density p(r) for / < L < M is given by ^ 

p{r) = -^^^ + 0{s + s = ^. (4.18) 

2[s + 1) rVl — T 71"/ 

The density has a cut-off for exponentially small r, which is irrelevant for {R±) if 7^ 1. 
Substitution of Eq. ( [4.18| ) into Eq. ( ^.17| ) yields the average reflectances 

= 2TRe °^ ° ^/ artanh a^, (4.19a) 
a*(a'^ - I) 

{R+) = 1- 2TRe ° ^/ artanh (4.19b) 

where T is again the transmission probability ( ^.10] ) from radiative transfer theory. Both 
quantities have a smooth L-dependence, with {R-) decreasing monotonically oc 1/L. In 
contrast, radiative transfer theory predicts a non-monotonic L-dependence for A > 1, leading 
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FIG. 6. Average reflectances {R±) as function of L/l for a = 7r/4 and 6 = 0.6, 0.9, in the 



coherent regime. The full curves are the analytical results for N ^ 1, computed from Eq. ( A18 ). 
The dotted curves are the large L/l limit given by Eqs. ( p^ and ( |42ll ). Data points are results 
from numerical simulations. 



to a divergence at some L. For A < 1, radiative transfer theory predicts a quadratic decrease 
of (-R-) oc l/L"^, for large L. The conclusion is that, in the coherent regime, phase coherence 
modifies the reflectance of the phase-conjugating mirror to leading order in N. 



The result ( 4.19 ) was obtained for the simplified model of a scalar reflection matrix 
a. The true a in Eq. ( [4 .61) is diagonal but not a scalar. This complicates the calculation 
because f2 = U]LaU+ then needs to be averaged over IA{N) even though U_ = U+. The 
calculation is outlined in App. 0. The complete result is a complicated function oi L/l 
(plotted in Fig. ^). For L/l ^ 1 the result takes the form of Eq. (|4.19| ), where now is to 
be determined from the equation 



^Tr ^ = 

1 — flO cL 1 — 



Tr = (4-20) 



In the limit N ^ oo this becomes an integral equation for qq, 

f"^/"^ COS0 a(0) ao 



ao a{(f)) I — Oq 



2 ' 



(4.21) 



where a{(j)) is given by Eq. (|2.6| ). As shown in Fig. ^ the large-L asymptote (|4.19| ), (f4.20|) 



is close to the complete result for L ^ I. In the limit 5 — *■ the solution to Eq. ( |4.21| ) is 
given by ao = 1.284 — 0.0133 z, for a = tt/4. The corresponding reflectances (for L ^ I) are 
= 61.1 l/L, and {R+) = 1 + 57.7 //L. 
To make contact with the work on wave-front reconstruction ||TT],|12|, we consider also 
the case of plane wave — rather than diffusive — illumination. A plane wave incident at 
frequency in mode n is reflected into modes m = 1, 2, . . . at frequency u}± with proba- 
bility (|(r±+)mnP)- The calculation of this probability proceeds similarly as the calculation 
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of R-. (See Ref. for the analogous calculation in the case of Andreev reflection.) Using 
Eqs. (|49|)- (|4.12| ) we can write 



r_+ = ie'^V*_OYl, (4.22a) 



O = ^T_0(1 - ^1 - - T_n)-'^T+. (4.22b) 

For the coherent regime, we may again identify V+ = V_ = V. The integration over U{N) 
can be performed using |]2^ 

KkVmjKiVml = j^2_ I ^^ik^Ol + ^mnkl^ji) " j^ i^kl^ji + ^mnkk^jl)- (4.23) 



We then find 

N + 1 m-N 



dy\{r-.Un\' = '4^R- + ^^I^EO.O*,. (4.24) 



In the limit of large we can write J2i^j OiiO*j = |Tr Op. In the same way as before, for 
L ^ I, this trace can be expressed in terms of oq, where Oq can be found form Eq. ([4.61): 
N~^Tt O = Tartanh oq. The result for the averages is then 

(|(r_+)„„p) = T^|artanhaoP, (4.25a) 
(|(r_+)^„H = N-'{R.), m^n. (4.25b) 

The incident plane wave is reconstructed with an intensity oc T^, in agreement with 
Refs. |Tl|JT^. In the coherent regime, off-diagonal (m ^ n) and diagonal (m = n) reflection 
probabilities differ by a large factor of order NT . 

In the incoherent regime, the matrices V+ and V_ are independent. Integration over 
IA[N) results in integrals of the form / dV V^fy*^ = N~^6ik- Then the off-diagonal and 
diagonal reflection probabilities are both given by 

(|(r-+)™„r) = Ar-i(i?_), (4.26) 

so there is no peak in the reflected intensity at the angle of incidence. This holds for every 
and L. 

For both the incoherent and the coherent regime we flnd for the reflection without fre- 
quency shift {uj^ —>■ uj^) the probability 

(l(r++)™nP) = \^(i?-,). (4.27) 

Here we see a much smaller backscattering peak, where the diagonal reflection probability is 
only twice as large as the off-diagonal reflection probability p6[ . This factor is independent 



of the phase-conjugating mirror, and exists entirely because of time-reversal symmetry [27 
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FIG. 7. Average reflectances {R±) as function of L/l for a = 7r/4 and 6 = 0.6, 0.9, in the 
incoherent regime. The curves are the analytical results for ^ 1, computed from Eq. ( ^.121) . 
Data points are results from numerical simulations. (Statistical error bars are shown when they 
are larger than the size of the marker.) 



C. Numerical simulations 



To test the analytical predictions of random-matrix theory we have carried out numerical 
simulations. The Helmholtz equation, 

(-V^ - eujl/c'^)S = 0, (4.28) 

is discretized on a square lattice (lattice constant d, length L, width W). Disorder is 
introduced by letting the relative dielectric constant e fluctuate from site to site between 



l±fc. Using the method of recursive Green functions we compute the scattering matrix 
Sdisorder(i^) of the disordered medium at frequencies uj+ and The reflection matrix rpcM 
of the phase-conjugating mirror is calculated by discretizing Eq. (|2.2|). From Sdisorder('^±) and 



rpcM we obtain the reflection matrix r of the entire system, and hence the reflectances ( [4.101) . 



We took W = 51d, Se = 0.5, a = ir/A, and varied 6 and L. For the coherent case we 
took = u;_ = 1.252 c/rf, and for the incoherent case = 1.252 c/d, uj_ = 1.166 c/d. 
These parameters correspond to = 22, /_|_ = 15.5 d at frequency u^. The mean free path 
is determined using Eq. ( ^.1U| ), which holds up to small corrections of order N~^. In the 
incoherent case we have = 20, Z_ = 20.1 d. This leads to Auj = 0.043 c/d and a dwell 
time for L/l ~ 3 of Tdweii — L/d ~ 150d/c. Hence we have TdwoiiAct; ~ 6.5, which should 
be well in the incoherent regime. For comparison with random-matrix theory, we take the 
large A^-limit and use the value for /. 

The numerical results are shown in Figs. |^ (coherent regime) and |^ (incoherent regime), 
for 6 = 0.6 and 0.9. As we can see the agreement with the analytical theory is quite 
satisfactory. The rapid rise of {R±) in the incoherent regime for the smallest 6 is accompanied 
by large statistical fluctuations, which make an accurate comparison more difficult. Still, the 
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FIG. 8. Histo grains for the niodcLl distribution |(r — [-)mi 
frequency shift for normal incidence. The results are for a single realization of the disorder at 
W = L = 251 d [L/l = 16.2), a = 7r/4, and 6 = 0.9. The arrow indicates the theoretical value 

(l(r~+)ii|' 



from Eq. ( [4.25| ), representing the ensemble average in the large A^-limit. 



striking differences between the coherent and incoherent regimes predicted by the random- 
matrix theory are confirmed by the simulations. 

We have also studied the backscattering peak for plane- wave illumination. We considered 
a square sample {W = L = 251 c?) with a = 7r/4, 6 = 0.9. We calculated the reflection 



probabilities 



for normal incidence {n = 1) in both the coherent and the incoherent 



regimes. The numerical results for a single realization of the disorder are shown in Fig. 



The arrow denotes the analytical ensemble average ( |4.25| ) of the backscattering peak in 
the large-iV limit, which is consistent with the numerical data. Notice the absence of a 
backscattering peak in the incoherent regime. 



V. COMPARISON WITH ANDREEV REFLECTION 

We have studied the reflection of light by a disordered dielectric medium in front of a 

The electronic 



phase-conjugating mirror. This problem has an electronic analogue [17,18 



disordered system consists of a metal, in which electron or hole excitations are scattered 
elastically by randomly placed impurities. Retro-reflection at the phase-conjugating mirror 
is analogous to Andreev reflection at the interface with a superconductor. The Fermi energy 
Ep plays the role of the pump frequency Uq, while the excitation energy E corresponds to 
the frequency shift Au. In spite of these similarities, the optical effects found in this paper 
have no electronic analogue. It is instructive to see where the analogy breaks down. 



To this end we compare the wave equation (|2.2| ) with the Bogoliubov-de Gennes equation 



H 

A* 



A 

-H 



E 



(5.1) 
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which determines the electron and hole wavef unctions u and v. The Hamiltonian 

H = -—V^ + V-Ef (5.2) 
Zm 

contains the electrostatic potential ^(r), which plays the role of the dielectric constant. 
(More precisely, ^0(5 — 1) corresponds to —{2m/lx^)V.) The role of the nonlinear suscepti- 
bility is played here by the pair potential A(r), which is only nonzero in the superconductor, 
where it equals Aqc"*'^. Comparing Eqs. ( p.lD and ( p.2[ ) for the electronic case with the 
optical equations ( |2.2|) and (|2.4|) one notices many similarities, and some differences which 
amount to a redefinition of quantities. There is however one essential difference: the matrix 
operator in Eq. ( ^.1| ) is Hermitian, while that in Eq. ( |2.2| ) is not, because of an extra mi- 
nus sign in one of the off-diagonal elements. This minus sign is the origin of the difference 
between Andreev reflection and optical phase conjugation. 

In the optical case a disordered medium becomes transparent (i?_ = 1) |9,10] for unit 



reflectance at the phase-conjugating mirror (a = 1). This does not happen in the electronic 
case, where R- is reduced by disorder even for ideal Andreev reflection. The reflection matrix 
of the normal-metal-superconductor (NS) interface, obtained from Eq. ( p.l| ) for = 0, 
E < Ao < Ep, is given by 

''NS = . n • (5.3) 



-le'^ 

Comparison with Eq. ([4. 61) for rpcM shows that Andreev reflection is independent of the 



angle of incidence; the matrix a in Eq. ( |4.6| ) is replaced by the unit matrix in Eq. ( [5.3| ). This 
is a substantial simplification of the electronic problem, compared with the optical analogue. 
The matrix tns is unitary, in contrast to rpcM, so that the appearance of gain or loss at 
the phase-conjugating mirror has no electronic counterpart. The reflectance R- = 1 — i?+ 
is a monotonically decreasing function oi L/l in the electronic case both in the coherent 
regime, 

i?„ = (2 + 4L/7r/)-\ if E < /i/Tdwcii and L > /, (5.4) 
and in the incoherent regime, 

= (1 + AL/tiI)-\ if E > /l/Tdwell- (5.5) 

The result ( |5.5| ) is what one obtains from Eq. ( |3.12| ) for the case A = 1 oi unit reflectance 
at the interface. (The transmittance T = (1 + 2L/ttI)~'^ of the disordered medium is the 
same for electrons and photons.) The result ( p.4| ) however, is not what one would expect 
from the optical analogue. Indeed, Eq. ( [4.17|) with ao = 1 would give R- = 1 for all L in the 



case of unit reflectance at the phase-conjugating mirror. The reason that the analogy with 
Andreev reflection breaks down is the difference of a minus sign in the wave equations ( p.2|) 
and (|5.1| ), which reappears in the reflection matrices ([4.6|) and (|573|) for phase conjugation, 



and ultimately in the reflectances in the coherent regime: 

R- = -— Tr ( 7 1 7^ 1 for electrons, (5.6a) 




i?_ = -— Tr =1 for photons if a = 1. (5.6b) 
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Here t and r are the transmission and reflection matrices of tlie disordered medium, wliicli 
satisfy tt''' + rr''" = 1. 

In conclusion, we have shown that the presence of a phase-conjugating mirror behind 
a random medium drastically changes the total reflected intensity, even when the medium 
is so disordered that wave-front reconstruction is ineffective. On increasing the frequency 
difference Acu between the incident radiation and the pump beams, a minimum in the 
disorder dependence of the reflected intensity appears. In a certain parameter range, the 
disordered medium reflects less radiation on reducing Atu. Experimental observation of this 
"darkening" would be a striking demonstration of phase-shift cancellations in a random 
medium. 
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APPENDIX A: CALCULATION OF THE REFLECTANCES IN THE 

COHERENT REGIME 



In Sec. 1^ we computed the average reflectances {R±) for the incoherent regime. For 
the coherent regime we presented only a derivation for scalar reflection matrix a. This 
appendix contains the calculation of {R±) for arbitrary matrix a. Our calculation is based 



on the diagrammatic method for integration over the unitary group of Refs. |^ and 
Integrals over the unitary group are needed for the computation of (-R±) because of the polar 
decomposition ([4.11|) of the transmission and reflection matrices. We find it convenient 
to use a slightly modified version of the diagrammatic technique, in which we apply the 
diagrammatic rules without making explicit use of the polar decomposition. We first outline 
the calculation of {R±) in which the diagrammatic method is used for the integration of 
the matrices U and V in Eq. ( [4.11|) , and then discuss the modification of the diagrammatic 



method. 

We start the calculation of {R±) with the elimination of the reflection matrix rii(ti;o) and 
the transmission matrices ti2(ti;o) and t2i(a;o) from the reflectances i?+ and R- [cf. Eqs. ( [4.9|) 



and ( |4.1CI| )], in favour of the matrix r = r22(wo). The result is 



{R+) = ^Tr [s'(a, a) - s(l, a) - s(a, 1) + s'(l, 1) + h(l) + h(l)*], (Ala) 

= i^Tr {s(a, a) - s'(l, a) - s'(a, 1) + s(l, 1) + - h(l) - h(l)*]}, (Alb) 

where we defined 

s'(x, y) = (x(l - rar*a)"^ryy*r*(l - a*ra*r*)"^x*), (A2a) 

s(x, y) = (x(l - rar*a)-^a-Vy*a"^*(l - a*ra*r*)-^x*), (A2b) 

h(x) = (x(l - rar*a)"^). (A2c) 



To perform the average over r, one may use the polar decomposition [cf. Eq. ([4.11|) ] 



20 



r = iVVT^V^, (A3) 

where U is a unitary matrix and T is the diagonal matrix containing the N transmission 
eigenvalues tj on the diagonal. The matrix U is a member of the circular unitary ensemble 
(CUE), i.e. it is uniformly distributed in the unitary group. The transmission eigenvalues 
Tj have density |^ 

p(r) = {2N/n) Imf/(l/r - 1 - iO,s), (A4a) 
f/(C,s) = cotanh[C-sf/(C,s)], s = 2L/ttI. (A4b) 

To integrate the matrix U over the unitary group, the matrices s, s', and h are first 
expanded as a power series in U. The integration of U is then done using the general 
expression for the average of a polynomial function of U [pO|] , 

n 

'Uaibi ■ ■ ■ Ua„6,„UQ,^^^ . . . Uq,^^^^^ = 6m,n X! ^l,---,Cfc Jl ^aj ,ap^j^Sbj ,l3p, ^^y (A5) 

P,P' j=l 

Here the summation is over all permutations P and P' of the numbers 1, . . . , n. The numbers 
Ci, . . . ,Cfc denote the cycle structure of the permutation P~^P'. (The permutation P~^P' 
can be uniquely written as a product of disjoint cyclic permutations of lengths Ci, . . . ,Cfc, 
with n = I]y=iCfc.) To compute {R±) in the limit of large N, it is sufficient to know the 



coefficients Vci....,cfe to leading order in A^. These are given in Refs. [^| and 1^^, together 
with a diagrammatic method which enables one to restrict the summation over P and P' to 
those permutations P and P' of which the contribution to {R±) is of maximal order in A^. 

Although the computation of {R±) is straightforward now, the actual calculation is rather 
cumbersome. We find it convenient to modify the approach of Refs. and ||2^ so that 



it can be applied directly to the average over the matrix r, without making explicit use of 
the polar decomposition (|A3|) . This is possible because the general structure (|A5|) already 



follows from the invariance of the distribution of U under transformations 

U ^ VUV, (A6) 

where V is an arbitrary unitary matrix. The fact that U itself is unitary is necessary to 
compute the value of the coefficients but it is not relevant for the general structure 

(|A5|). Since the matrix r is both unitary and symmetric, its distribution is invariant under 
transformations 

r ^ VrV"^ (AT) 

that respect the symmetry of r. The same group of transformations leaves invariant the 
circular orthogonal ensemble (COE), consisting of uniformly distributed unitary and sym- 
metric matrices. A diagrammatic technique for averages over the COE is presented in Ref. 



24 1 . As before, the general structure of the average of a polynomial of a matrix from the 



COE is entirely determined by the invariance under the transformations (|A7|) , and therefore 



applies to the reflection matrix r as well. It reads p^pl 



2n 
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" trh(z) * trh(z) tr h (z) 


FIG. 9. Diagrams for the calculation of ho(z). 



where now the summation is over permutations P of the numbers 1, 

P as 



2n. We may write 



(A9) 



where the permutations aj and a'j operate on the numbers 2j — 1 and 2j, and the permutation 
Pe (Po) permutes even (odd) numbers only. The numbers Ci, . . . , in Eq. ( |A8| ) are the cycle 
structure of the permutation P~^Po. The specific values of the coefficients Vci,,,,,Cfc for an 
average of r are of course different from those for the COE. 

Now that we have identified the formal equivalence of an average over the (non-unitary) 
symmetric reflection matrix r and a unitary symmetric matrix from the COE, we can directly 
apply the diagrammatic rules of Refs. ||29| and p4| to an average over the matrix r, provided 



we know the coefficients V^i,...,cfc for the ensemble of reflection matrices r of a disordered 
waveguide. To flnd these coefficients, we use the fact that they factorize, to leading order 
in N, 



Cl,...,Cfc 



just as they do for the COE. This follows directly from the fact that 
in A^, the average {UjTr{rr^p) factorizes into Uj{Tr (rr^Y^) H 
coefficients K. Hereto we consider the function 



(AlO) 

to leading order 
It remains to find the 



ho(^) 



1 — rr*z'^ 



(All) 



We first compute h^lz) from the diagrammatic technique, with a priori unknown coefficients 
Vc- We then compare our result with a calculation of Tr ho(-2) from the density of transmis- 
sion eigenvalues ( |A4]) . The relevant diagrams for the diagrammatic calculation are shown 
in Fig. ^ (for a detailed explanation of the diagrammatic notation of Fig. ^ we refer to 
Ref. ||2^). The result is a self-consistency equation for ho(-2) that involves the generating 
function F of the coefficients K, 
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' l-zF[TTho{z)y ^^^^^ 

oo 

Fix) = J2V,x^'-\ (A13) 

j=c 

Here IL is the N x N unit matrix. Direct computation of Trho(z) from the density p(r) of 
transmission eigenvalues, gives 



1 p{t)z 
' I-z^(I-t) 



Trho(;.) = / dr ^ f _^ . (A14) 



Together, Eqs. ( |A1U| )-( [A14D determine the coefficients K;i,...,cs, needed for the diagrammatic 



evaluation of {R±)- In the limit of L — oo, the density of transmission eigenvalues tends 
to N6{t). Hence Iiq^z) = z/{1 — z"^) and F{x) = {y/N"^ + Ax"^ — N)/2x. The corresponding 
coefficients Vc = N^^"^^ {^^^f^ are precisely those of the COE For finite L, the density 
p(r) is no longer a delta function, and hence the coefficients Vc deviate from those of the 
COE. 

The fact that we can use the diagrammatic rules directly for the average over r simplifies 
the calculation considerably. A central role in the calculation is played by the function h(x) 
defined in Eq. (|A2|) . The diagrams for the calculation of h(x) are similar to those of Fig. 
and the result is a self-consistency equation for h(x) 

h(x) = x(l - aF[Tr h(a)]) "\ (A15) 



Notice the formal equivalence with Eq. ( |A12[ ). The function F was defined in Eq. ( [A IS] ) 



Using the diagrammatic technique for the computation of s and s', we find the linear relations 



s(x,y) = h(x) 



— 1 * — 1* 

a yy a 



K Tr s'(a, y) + L Tr s(a, y)! h(x) 



s'(x, y) = h(x) [K Tr s(a, y) + L Tr s'(a, y)] h(x)*, 
where we have defined 



/i = Tr h(a) 

oo 



^ = 1^ Vi+j-ih [h = _ 



h*F{h) - hF{h)* 

J^2 _ J^*2 



(Al6a) 
(Al6b) 



(Al7a) 
(A17b) 

(A17c) 



The relevant diagrams leading to Eqs. ( |A16| ) and ( |A17| ) are shown in Fig. |10[ They are 



similar to those of Ref. where the case of a chaotic cavity was considered, instead 
of a disordered waveguide. Together, Eqs. ( |A15| )-( |Al7[ ) form a closed set of equations, 
from which s(x, y), s'(x, y), and h(x) can be calculated. The average refiectances {R±) are 
obtained upon substitution of s(x, y), s'(x, y), and h(x) into Eq. (|AT|). The final result is 
expressed function of h = Trh(a), 
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FIG. 10. Diagrams for the calculation of s(x, y) and s'(x, y). 
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N{R^) ^ + -^''/^^, ^'[gj)^"' + ^R^Tr [1 - af (A)]-', (A18a) 
where we defined 

/ = Tra(l - aF{h))-\l - a* F{hy])-^a* , (Al9a) 

J = Tr (1 - aF{h))-\l - a*F{h))-\ (Al9b) 

These expressions simplify in the large L//-hmit, when p(r) takes the form ([4.18|) . Sub- 
stitution in Eq. (|A14|) gives 

, / N Nz ( zartanhz\ ,n ^ , , /. n 

TrhoZ=^ 5 1 ^ 1 + s = 2Lnl. A20 

and hence allows us to find F{z) from Eq. (|A12|) . Expanding the expressions (|A18|) for {R±) 
and the self-consistency equation ( [A15| ) to lowest order in (l + s)"^ we find the results ( [4 .191) 
and (|4.2CI| ), with the effective reflectance Oq = z. 
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